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Abstract. Recent numerical and observational studies revealed that spectra of mag- 
netic and velocity fluctuations in MHD turbulence have different scaling indexes. This 
intriguing feature has been recently explained in the case of weak MHD turbulence, that 
is, turbulence consisting of weakly interacting Alfven waves. However, astrophysical 
turbulence is strong in majority of cases. In the present work, we propose a unifying 
picture that allows one to address weak and strong MHD turbulence on the same foot- 
ing. We argue that magnetic and kinetic energies are different in both weak and strong 
MHD turbulence. Their difference, the so-called residual energy, is spontaneously gen- 
erated by turbulence, it has the Fourier spectrum E r (k) = E v (k)-Et,(k) oc -f w (k\\/k ± )k^ 2 
in weak turbulence, and E r (k) oc -f s (k\\/k ± )k2 3 in strong turbulence. Here f WiS (x) are 
functions declining fast for x > C w>s and not significantly varying for x < C W)J with 
some constants C w> „ and k\\ and k ± the field-parallel and field-perpendicular wave vec- 
tors with respect to the applied strong uniform magnetic field. 



1. Introduction 

Magnetic fields and turbulence are common in a variety of astrophysical plasmas, from 
planets and stars to interstellar and intergalactic media. Magnetic turbulence is also 
commonly invoked to explain small-scale features of the solar wind. Numerical simu- 
lations and analytic modeling play an important role in interpreting observational data. 
Recently, it has been found that magnetic and velocity fluctuations are no t in equipar- 
tition i n MHD turbulence ( e.g., iPodesta et ail 120071 : iTessein et all 120091 : IChen et all 



l2011al : lBoldvrev et all 1201 lh . which seems to be at odds with basic assumptions of 



conventional models of MHD turbulence. The goal of the present contribution is to 
propose an explanation for this intriguing phenomenon. In contrast with ordinary in- 
compressible turbulence, which is always in a strongly coupled state, incompressible 
MHD turbulence can exhibit two distinct regimes of weak and strong turbulence. This 
stems from the fact that the MHD system possesses Alfven waves that can coalesce 
and scatter off each other. When during a singe interaction the wave amplitudes change 
only slightly, turbulence is weak, otherwise, it is strong. It is important to note however 
that even if MHD turbulence is weak at large scales, its strength increases toward small 
scales, so that the range of scales where weak MHD turbulence may exist is typically 
limited. In this contribution we present a unifying model of MHD turbulence valid 
for both weak and strong regimes. The incompressible MHD equations for magnetic 
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and velocity fields, b(x, t) and v(x, t), have especially useful form when written in the 
so-called Elsasser variables z ± = v + b, 

j t +y A -Vjz ± + (z r -V)z ± = -VP. (1) 

The equations are written in a frame with zero mean-flow velocity, b is the fluctuating 
magnetic field normalized by ^j4npo, va = Bq/ -\jAnpo is the Alfven velocity corre- 
sponding to the uniform magnetic field Bo (so that the total magnetic field is B = Bo+b), 
P = plpo + b 12, it includes the plasma pressure, p, and the magnetic press ure, pq is 
the constant mass density, and we neglect driving and dissipation terms (e.g., Biskamp 
20031 : lMarsch~& Mangeney 1987). In what follows we assume that turbulence is driven 



at large scales, that can be mimicked by adding forcing terms to the right-hand sides of 
Eqs. (Q]). Small-scale turbul ence is expected to be independent of the large-scale driv- 
ing (e.g. jMason et al.ll2008h . We will also assume that the uniform guide field is strong 
compared to the rms fluctuations, that is, b rms «; Bq. 



2. The energy spectrum 

The ideal MHD equations conserve the two Elsasser energies, E ± = J^^d^x - 

f e ± (k)d 3 k. When the energies are supplied to the system at large scales, they get re- 
distributed over scales by nonlinear interactions, and removed from the system at small 
dissipation scales. One can argue that the energy gets redistributed predominantly over 
the modes whose wevevectors are approximately normal to the strong guide magnetic 
field. We will concentrate on the so-called balanced case, when e + ~ e~, and we can 
therefore represent the Fourier energy spectra in the form 

e ± (*||,*J = / ± (*||/*x)*r ,r . ( 2 ) 

where / ± (jc) do not vary significantly for x < C and decline fast for x > C, with some 
constant C. Here k\\ is the wavevector in the direction of the uniform field Bo, and k x is 
the wavevector in the field-perpendicular direction. This form of the spectral functions 
is motivated by the fact that the energy redistribution occurring due to small-scale fluc- 
tuations (large k) is predominantly normal to the direction of the local guide field, which 
is the field produced by large-scale fluctuations. Therefore, compared to the direc- 
tion of the global uniform field, the energy spectrum is smeared inside the small angle 
#o ~ b rms /Bo, which implies a wedge-shap ed energy-containing domain k\\ < Oqk , , or. 



the spectral functi on (O with C ~ Gq, (e.g. JCho & Vishniaq|2000l : iMaron & Goldreichl 



l200ll : IChen etafll201 lbh . One can then write down a model equation for the spectral 



function ©, using certain closure assumptions. Differentiating e ± (k) with respect to 
time, iterating equation CD once, and splitting the forth-order correlation functions of 
z's into pair-wise correlations using Gaussian rule, we gefl 

dtfQm, *J - J M(k, ki, k 2 )0 ± (/c 2 ||, k 2± )e T {k 2l] , k 2± ) [^(tin.^ij.) - e ± {k h k A _)} x 



'Some extra assumptions are made in obtaining this equation, for instance it is assumed that the cross- 
correlation (z + ■ z~) is absent, see, e.g., Goldreich & Sridhar d 19951) . It should however be borne in mind 
that this equation is not rigorously derived from l[T); it should be considered only as a model equation or 
as a plausible two-point closure. 
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x6{k\\ - ky\\ - k2\\)S(k± - ki ± - k2 ± ) d 3 k\ d 3 ki- (3) 

In this equation, the kernel has the form Mk,k,,k 2 x (k± x k2±) 2 (k ± ■ k\±) 2 /(k?[_k% ), 
and the 0* functions depend on the assumptions about the nonlinear interaction made 
in the model. In general we argue that these functions should concentrate in the region 
where the nonlinear interaction is essential, and inside this region they should scale as 
the inverse time of nonlinear interaction, that is, ± (^||,^ ± ) oc l/r(k±). This can be 
summarized as follows, & ± (k\\,k ± ) = g ± (k\\,k A _)k~ i s , where g ± (k\\,k±) « const in the 
region of nonlinear interaction, and the nonlinear interaction time scales as r(k ± ) oc k s ± . 

To understand better our model ([3]>, consider particular examples. In the case of 
weak turbulence, g(k\\,k ± ) « const in a quite narrow region k\\VA < 1/t(& ± ) compared 
with the ^||-widths of the functions e ± , and it declines fast for 1c\\Va > l/f(k±). There- 
fore, &(k\\, k ± ) is a broadene d delta-function of h , obeying J &(k\\, k±)dk\\ = const. One 
then recovers the theory by iGaltier et all (l2000h . In the case of strong turbulence, one 



expects the nonlinear interaction to be important in the same region where the energies 
are concentrated, that is, g* « const in the region k\\ < 8ok t , it declines outside 



of this region, and r ~ A/z(sl). This way we recover the iGoldreich & S ridhar ( 19951) 



theory. If in addition to the assumptions of the GS theory, one assumes that there is per- 
sistent dynamic angular alignment between magnetic and velocity fluctuations, which 
reduces the nonlinear interaction by Ox ~ /t 1 ^ 4 , one needs to multiply the kernel in © 
by 9\ ~ k~ 1/2 , and assume that r ~ A/{vx6a)- O ne then recovers the theory by lBoldyrevI 
(200% In view of this, we stress that model © provides a useful description of the 
spectral energies in MHD turbulence, however, it crucially depends on the scaling as- 
sumptions about the interaction time, incorporated in the modelQ 

The steady spectrum of turbulence can then be found by requiring that the col- 
lision integral in the rhs of j3]) i s zero. This leads to the spectrum of weak turbu- 
lence (|Ngt&lhattach^e3 Qjwl; |GaltieretaD|2QQ3): E ± (k h k ± ) = e ± {k\\,k A _)2nk 1 _ oc 
f*(k\\)k~ 2 , where f^(k\\) depend o n the details of large-scale dri ving. The field-perpendicular 



spectrum of strong turbulence in Goldreich & Sridhari (Il995h theory is then found as 



E ± (k ± ) = J e ± {k h k A _)2nk^dk\\ oc k ± 5/3 , while the spectrum in iBoldvrevI (l2006h the- 
ory is E ± (k ± ) = J e ± (k\\, k ± )2nk ± dk\\ oc k" 3 ^ 1 . Numerical simulations do produce the 
spectrum k~ 2 for weak turbulence ( Perez & Boldvrevll2008h. and the spectrum k~ 3 ^ 2 for 



stron g turbulence with a strong guide field Bo, e.g., (IMtiller & Grapp in 2005; M ason et al.l 



2008) 



3. The spectrum of the residual energy 

Recently, it has been realized that significant role in turbulence dynamics is played 
by the so-called residual energy, that is, t he energy difference between magnetic and 
kinetic fluctuations, E r = E v - Et,, see (Boldyrev & Perez 120091 : IWangetaD 120 111 : 
iBoldyrev et al.ll20"iT ). Indeed, the complete description of the second-order statistics of 
two fluctuating fields, v and b requires three independent correlation functions. Two of 
them are provided by the autocorrelation functions of the Elsasser variables, that is, the 



2 The same statement is true for the so-called EDQNM closures often used to derive equations of type (TJJ 
for the spectra of strong turbulence. While providing physically reasonable models of turbulence, such 
equations are not derived from first principles and they crucially depend of model assumptions. 
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energy spectra ©. The third one, the cross-correlation function, is the residual energy 
E r = J (z + • iT) (fix = J Re[e r (k)] d 3 k, where, by definition, e r (k) = z + (k) • z"*(k). It 
is easy to see that e' (k) is a complex function, while the residual energy spectrum is its 
real part. 



1976; 



2007: 



The residual energy has been previously addressed in the literature (e.g., Pouquet et al. 
Grappin et al.ll983l : IZank et all 19961 : iMuller & Grappinll2005l : lNg & Bhattacharied 
Chen et al. «201 lah . ut it has been studied to a much lesser extent compared to the 



Elsasser energies, possibly because it is not a conserved quantity, it is not sign-definite, 
and and it does not exhibit a cascade in a turbulent regime. As a result, it cannot be 
expressed through the conserved Elsasser energies (or, equivalently, through the total 
energy or cross-helicity), which are commonly used in theoretical models and mea- 
sured in observations. In fact, in many studies of MH D turbulence the res idual energy 
is explicitly or implicitly assumed to be zero, see e.g., ( Galtier et al.ll200Ci|) . In this sec- 
tion we propose a model for the residual energy, analogous to Eq. Q. We demonstrate 
that in contrast with the scaling of the Elsasser energies, the scaling of the residual en- 
ergy is quite robust, that is, it depends to a mush lesser extent on the arbitrary scaling 
assumptions incorporated in the model. 

To obtain the equation for the residual energy, we first note that in the absence 
of the nonlinear interaction, the spectral residual energy e r (k) = z + (k) • z~*(k) oscil- 
lates in time, since z + (k) oc exp(ik\\VAt) and z~(k) oc exp(-/&||VA?)- When the nonlinear 
interaction is present, the residual-energy evolution equation should contain the terms 
describing interaction of the residual energy with the Elsasser fields, - e r e ± , and gener- 
ation of the residual energy by the Elsasser fields, ~ e^e~ . It has been recently realized 
that the terms descri bing the generatio n of the residual energy by the Elsasser fields are 
essentially nonzero (IWang et al.ll201 ll) . meaning that residual energy is spontaneously 
generated by turbulent dynamics even if it is zero initially. We start our discussion of 
the residual energy with more detailed consideration of these terms. 

It is crucial to note that the terms describing generation of the residual energy by 
the Elsasser fields should have the same dimension as the rhs of Eq. Q. Indeed, the 
residual energy has the same dimension as the Elsasser energies, and it is generated due 
to same nonlinear interactions. We however do not need the exact structure of this term, 
rather, we need to know its scaling with respect to the wavenumber. It turns out that this 
scaling is rather universal. Indeed, the term in the rhs of Eq. (f3]> describes constant flux 
J of the Elsasser energies e ± (k ± ) in the field-perpendicular direction, that is, it should 
scale as x ^j 2 * no matter what particular model of turbulence is assumed. The 

term describing generation of the residual energy should then have the same scaling, 
although its structure is different. We therefore model the residual-energy generating 
term as a(k\\,k ± ) with the only requirement that it is concentrated in the region where 
the Elsasser energies are concentrated, and it obeys f a(k\\,k±)dk\\ oc &~ 2 . 

The term describing relaxation of the residual energy due to its interaction with the 
Elsasser energies can be generally modeled as —y(k\\, k ± )e r (k\\,k ± ), where the relaxation 
rate y depends on the spectrum of the Elsasser fields, and it is concentrated in the 
region where the nonlinear interaction is present. We now collect all the three terms to 
formulate our model equation for the residual energy: 



d t e r (k\\,kj_) = 2ik\\VAe r - y{k\\,k A _)e r + a(k\\,k±). 



(4) 
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4. Discussion 



We now apply our formalism to the cases of weak and strong turbulence. For weak 
turbulence, the energy cascades predominantly in the field-perpendicular direction for 
each k\\, so that the field-parallel structure of the spectrum does not change with k±. 
Moreover, in this c ase the residual-energy generating term can be shown to be negative 
dWang et alJIioTlh . We can therefore write (restoring the dimensional coefficients) that 



a(k\\,k ± ) - -a w {v^ m JvA)k~] 1 for k\\ < ko\\, where a w is a dimensionless constant and 
&o|| ~ is the field-parallel spectral width of the Elsasser fields. Weak turbulence 
theory also allows one to estimate the time of nonlinear interaction of the fields, which 
gives y - fi(v 2 m JvA)k± for k\\ « (it will be clear momentarily why only the region 
k\\ ~ is relevant here), and ft is a dimensionless constant. Solving Eq. © for long 
times, we get: 

e v -e b = Re[e (k lh k x )] = (5) 

where e = v rms /vA ~ b rms /Bo «: 1. The residual energy is conce ntrated in a narrow 
region around kn < /3e 2 k i /2, in agreement with previous findings (Boldyr ev & Perez 



2009: [Wang et aljfeoi 1 ). The phase-volume compensated field-perpendicular spectrum 



of the residual energy then has the structure 



E r {k h k ± ) = Re[e r {k h k ± )]2nk ± - -M^/kJk- 2 , (6) 

where f w (x) = v 2 ms a w fie 4 /(fi 2 e 4 + Ax 2 ), as follows from (f5]). We can also define the 
field-perpendicular spectrum 



■(k\\,k ± ) dk\\ = a w nv r 



c 6 k , 



-V 2 € 2 k- 



(V) 



For the case of strong turbulence, the spectra of e ± (k\\, k ± ) are concentrated in the region 
k\\ < 0ok ± , and it is reasonable to assume that the function a(k\\,k±) is concentrated 
in the same region. Restoring the dimensional parameters, one can therefore write 
a(k\\,k±) = -a s (vl ms /L ± )(8ok ± )~ 1 k']^ for a given k\\ inside the region k\\ < 6ok ± , where 
a s is a dimensionless parameter, and L ± is the integral field-perpendicular scale of the 
fluctuations. Note that the power of k ± is fixed by the requirement J a(k\\, k ± ) dk\\ oc 
k~ 2 . One can also assume the power-law behavior for the relaxation rate, y{k\\,k ± ) = 
yk± within the same region, where y is a dimensional parameter. It should be noted, 
however, that in contrast with the function a(k\\,k ± ), whose scaling could be established 
on dimensional grounds, the scaling of the function y(k\\, k x ) cannot be easily derived. 
One can only argue that this relaxation rate should compete with the linear frequency 
only in the region where the Elsasser energies themselves are concentrated, that is, 
/i < 1. It is interesting, however, that this bound is enough to establish the field- 
perpendicular spectrum of the residual energy. The solution of (0]) takes the form 



(8) 
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To find the field-perpendicular energy spectrum, one integrates this result over k\\. One 
does not need however to integrate the function a(k\\,k±), since the prefactor is a nar- 
rower function (a broadened 5-function, in fact). The integral over k\\ is then indepen- 
dent of yW, and the result is: 

E r (k ± ) = J 'Re[e r (k h k ± )]2jrk ± dk\\~ -v 2 rms L~ ± l k' 2 , (9) 

where we used On ~ v r m JvA- This result is in agreement with numerical studies (e.g., 
Mull er & Grappinll2005h . It is also of interest to establish the value of //. This can 

be inferred from numerical simulations if one evaluates E r (k\\ = 0,k ± ) oc k 2 M . Our 
simulations (that will be reported elsewhere) indicate that, quite interestingly, ji « 1, 
which allows us to write the general expression for the residual energy in the form: 

E r (k h k ± ) = Re[e r {k h k x )Jlnk x = -f s (k\\/k x )k ± 3 , (10) 

where f s (x) * l/0n for x < On, and it declines for x > On. Relations @, ((6]|7]>, and d9]- 
[TOl are the main results of our work; they provide a model for residual energy observed 
in the solar wind and in numerical simulations. 
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